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1.  INTRODUCTION 


Researchers  involved  in  product  development  regularly  desire  to  modify  or  extend  computer 
models.  Usually  these  changes  or  applications  were  unforeseen  when  the  code  was  originally 
developed.  Altering  the  code  to  accommodate  the  researcher  represents  the  obvious  and  tradi¬ 
tional  solution.  For  a  large,  complex  code  this  represents  a  formidable,  and  frequently  frustrat¬ 
ing,  task.  The  situation  is  exacerbated  if  the  code  authors  are  unavailable,  the  algorithm  is 
poorly  documented,  the  code  is  unstructured,  etc.  Furthermore,  some  codes  are  extremely  fra¬ 
gile  due  to  obsolete  programming  practices.  Seemingly  small  or  insignificant  changes  often  will 
break  the  code. 

Proprietaiy  codes  represent  another  difficulty.  In  some  situations  source  code  is  specifically 
not  provided  for  a  variety  of  reasons.  This  scenario  usually  dictates  a  negotiated  contract  to 
implement  new  changes  or  features.  Besides  being  untimely,  the  process  leaves  little  margin  for 
error  or  dampens  the  desire  to  experiment  with  new  ideas. 

Faced  with  these  difficulties,  the  program  is  rewritten  more  often  than  not.  In  the  mean¬ 
time,  some  calculations  are  avoided  altogether  or  are  performed  via  tedious  cycles.  There  are 
alternatives  to  this  expensive  proposition.  It  is  quite  possible  to  utilize  an  established  computer 
model  as  an  isolated  black  box  and  write  another  code  which  merely  calls  it  as  a  programming 
sub-unit.  The  method  has  many  advantages.  Not  directly  tampering  with  the  innards  of  the 
black  box  is  just  one.  Another  strength  is  that  it  tends  to  modularize  the  problem.  Changes 
are  so  easy  to  understand  and  implement  that  it  encourages  further  investigation  and  algorithm 
refinement. 

This  report  will  illustrate  this  technique  using  the  popular  interior  ballistics  codes  IBHVG2 
[1]  and  BLAKE  [2].  Some  problems  will  be  posed  and  solved.  These  case  studies  were 
designed  to  be  instructive,  so  they  were  deliberately  simple  constructions.  Still,  we  feel  th^ 
represent  nontrivial  extensions  to  the  respective  codes  and  are  typical  of  calculations  which  are 
currently  performed  by  hand.  The  methodology  is  quite  general,  however.  The  reader  is 
encouraged  to  expand  on  these  examples  and  think  up  new  applications. 


n.  IMPLEMENTATION 


Calling  one  compiled  code  (already  executable)  from  another  can  usually  be  accomplished 
via  system  calls.  Since  the  Department  of  Defense  has  largely  converted  to  the  UNIX  operat¬ 
ing  system,  our  example  will  explicitly  cover  this  circumstance.  Other  operating  systems  will 
typically  allow  the  same  approach  but  the  details  will  vary. 

Figure  1  diagrams  the  process.  The  main  program  drives  the  problem.  In  our  simple  exam¬ 
ples,  it  declares  storage,  initializes  variables,  and  calls  a  subprogram,  SUBl,  to  perform  the 
mathematics.  This  program  makes  a  system  call  which  runs  MAIN2,  our  “black  box”  program, 
as  a  subordinate  unit.  The  results  of  this  run  are  available  to  SUBl.  We  hope  the  complete 
case  examples  which  follow  will  clarify  the  technique. 


1 


MAIN 


MAIN2 


**  Black  Box'* 
Program 


Execute  new 
process  via 
UNIX  shell 
script. 


Cominands 


Figure  1.  “Black  Box”  Flowsheet 
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A.  An  Optimization  Example 


A  charge  designer  frequently  must  maximize  the  projectile  velocity  given  a  set  of  con¬ 
straints.  Usually,  a  number  of  variables  are  available  which  may  be  manipulated  to  increase  the 
velocity.  Propellant  mix,  grain  dimensions,  and  charge  weight  are  typical  examples.  Velocities 
which  result  in  a  violation  of  a  design  constraint  are  infeasible.  As  an  illustration,  a  set  of 
inputs  resulting  in  a  breech  pressure  exceeding  safety  specificationswill  be  excluded.  This  appli¬ 
cation  is  a  classic  constrained  optimization  problem. 

The  example  problem  involves  optimizing  the  performance  of  a  120-mm  gun  system  using  a 
19-perf,  granular  propelling  charge.  For  this  system,  propellant  characteristics  are  varied  to 
maximize  the  velocity.  The  grain  is  characterized  by  setting  the  length  and  grain  diameter 
(assuming  perforation  diameter  to  be  fixed).  The  remaining  free  variable  is  the  total  charge 
weight.  Ail  other  system  parameters  remain  fixed.  The  constraints  are  easy  to  understand. 
Clearly,  the  designer  cannot  specify  more  propellant  than  can  actually  fit  in  the  gun  chamber. 
How  much  you  can  load  will  be  very  dependent  on  grain  dimensions,  however.  In  other  words, 
for  a  set  of  inputs  to  be  feasible,  the  charge  weight  must  not  exceed  the  maximum  loadable 
weight  for  a  given  grain  type.  Furthermore,  the  resulting  maximum  breech  pressure  must  fall 
below  the  given  design  specification. 


Table  1.  Notation  for  Optimization  Problem 


Symbol 

IBHVG2 

Variable 

Description 

1 

LEN 

grain  length 

D 

DIAM 

grain  diameter 

C 

CHWT 

propellant  charge  weight 

C* 

- 

maximum  loadable  charge  weight 

K 

CHAM 

chamber  volume  available  for  propelling  charge 

P’ 

- 

maximum  breech  pressure  for  a  given  set  of  inputs 

Pm 

PMAX 

maximum  allowable  breech  pressure  (design  variable) 

V 

VMUZ 

muzzle  velocity  of  projectile 

Using  the  notation  of  the  previous  table,  this  problem  can  be  simply  stated  mathematically 


as 


Maximize:  v  =g{l,D) 


subject  to  the  constraints 


Pu>P’  [l,D,C  I  C<C*(/,D)) 


(A.1) 


(A.2) 


Some  simple  sketches  should  clarify  the  constraints.  Maximum  breech  pressure  is  assumed 
to  be  a  monotonicaUy  increasing  function  of  charge  weight  for  a  given  grain  type.  In  this  case, 
there  will  be  two  possibilities.  Figure  2  illustrates  the  situation  for  a  given  grain  size  (/,D) 
where  the  maximum  loadable  charge  weight,  C’,  results  in  a  maximum  pressure  below  the 
design  limitation.  Figure  3  depicts  the  case  where  the  pressure  limit  dictates  a  loading  below 
C’ .  If  one  further  assumes  that  muzzle  velocity  is  a  monotonicaUy  increasing  function  of  max¬ 
imum  pressure  for  a  given  grain  type,  then  the  optimization  domain  would  not  include  C  since 
it  is  not  a  free  variable.  The  optimal  C  would  clearly  be  C*{l,D)  for  the  first  case  and 
f~^  (P^)  for  the  second  case  given  that  ?*=  / (C). 


Table  2.  Program  Organization 


Figure  1 

Example  1 

Description 

MAIN 

OPTIMZ 

main  program  for  optimization 

SUBl 

FUNC 

user-written  routine  caUed  by 
ZXMIN  (IMSL  routine) 

sys_ 

SYS 

C-code  to  perform  system  caU 

Commands 

run.ibhvg2 

operating  system  commands  to 

execute  “black  box” 

MAIN2 

IBHVG2 

“black  box”  code 
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c  (I,D) 


Figure  2.  Maximum  Loading  Limitation 
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Now  we  will  discuss  the  actual  coding.  In  our  specific  case,  we  desire  to  optimize  projectile 
velocity  subject  to  the  constraint  of  79  Kpsi  for  a  maximum  pressure.  The  manipulative  vari¬ 
ables  will  be  grain  diameter  and  length.  All  other  inputs  which  define  the  interior  ballistics 
remain  constant.  This  problem  was  solved  using  a  small,  user-written  program,  a  UNIX  shell 
script,  and  the  International  Mathematical  Subroutine  Library  (IMSL)  [3].  We  will  discuss  each 
in  detail 

The  main  program,  OPTIMZ,  primarily  does  one  thing:  it  calls  the  optimization  routine, 
ZXMIN,  from  IMSL.  It  is  quite  short  because  ZXMIN  performs  the  real  work.  This  was 
intentional.  IMSL  is  a  commercial  package  popular  to  many  scientific  organizations.  The  rou¬ 
tines  are  supported  by  a  professional  staff  and  are  extensively  tested.  We  chose  IMSL  because  it 
is  nearly  ubiquitous;  however,  there  are  many  other  high  quality  libraries  to  pick  from.  Some 
recommended  alternatives  are:  NAG,  MINOS,  MINPACK,  UNPACK,  CMLIB,  TOMS, 
Harwell  PORT,  and  SLATEC  (see  [4]  for  more  information).  Whatever  the  choice,  we  feel 
that  if  at  all  possible,  the  engineer  should  avail  himself  of  an  immense  amount  of  work  by 
merely  linking  to  pre-tested  routines.  Chosing  the  right  algorithm  for  the  job  can  be  challanging 
in  itself;  conducting  original  research  in  optimization  theory  is  usually  counter-productive  and 
better  left  in  the  hands  of  an  expert.  Furthermore,  the  documentation  and  verification  process 
is  simplified  because  one  may  reference  the  corresponding  library. 


OPTIMZ 


C 

c 

c 


c 

c 

c 


program  optimz 


external  func 

integer  icr,  n,  nsig,  iopt,  maxfn 
real  g(2),  h<13),  v,  x(2),  w(16) 

Initialize 


n  a  2 
nsig  a  3 
maxfn  a  lOO 
iopt  a  3 
x(1)  a  0.357 
x(2)  a  0.408 


#  dimension  of  domain 

#  number  of  significant  digits 

#  max  nunber  of  function  evaluations 

#  options  selector  (See  IHSL  manual) 

#  grain  diameter  —  initial  guess 

#  slot  size  —  initial  guess 


IMSL  routine  for  function  minimization 


call  zxmin(func,  n,  nsig,  maxfn,  iopt,  x,  h,  g,  v,  w,  ier) 


write  (6,^<f10.5,f10.5,f10.5)^)  x(1),  x<2),  v 
write  (6,'(a,i3)^)  Mer  a  » jer 
e^ 


ZXMIN  uses  a  quasi-Newton  method  to  find  the  minimum  of  scalar  function  of  N  variables. 
We  found  that  computing  the  Hessian  matrix  produces  a  conservative  but  controlled  optimiza¬ 
tion  (IOPT  *3)  for  this  type  of  problem.  It  should  be  added  that  there  are  numerous  optimiza- 


6 


tion  algorithms  available.  We  used  ZXMIN  merely  for  convenience.  No  effort  was  made  to 
customize  an  algorithm  for  the  type  of  topology  found  in  ballistic  problems.  Some  useful  alter¬ 
native  candidates  might  be  the  downhill  simplex  method  or  one  of  the  direction  set  methods 
[5,6,7].  These  should  execute  much  faster  because  no  derivatives  are  estimated  numerically. 

ZXMIN  requires  a  user-supplied  subroutine  to  compute  the  scalar  function  to  be  minim¬ 
ized.  FUNC  shown  below  performs  our  particular  task.  Since  we  desire  to  maximize  the  velo¬ 
city,  we  returned  the  negative  value  to  ZXMIN.  FTJNC  is  rather  unorthodox,  but  forms  the 
crux  of  this  report.  FUNC  serves  as  the  bridge  to  the  “black  box”  code,  IBHVG2.  The  role  of 
the  subroutine  is  to  generate  a  string  from  the  arguments  passed  into  it,  and  with  that  string, 
execute  a  UNIX  shell  script.  Our  Fortran  compiler  does  not  allow  such  a  system  call,  so  we 
link  in  a  C  program  called  SYS  which  will. 


FUNC 


SUBROUTINE  FUNC(N,  X,  V) 

c  call  a  UNIX  shell  script  which  executes  the  IBHVG2 
c  computer  code  with  two  parameters: 
c 

c  X(1)  —  grain  length 

c  X(2>  —  grain  diameter 

c 

c  and  returns  the  value  V  (muzzle  velocity) 


INTEGER  H 
REAL  X(N),V,P 
CHARACTER*! 20  CMO 
CHARACTER*25  ST0RE1 
CHARACTER*2S  ST0RE2 


c  using  internal  writes  to  convert  the  floating  point 
c  nunbers  to  a  character  string 


WRITE(  ST0RE1,  '(f25.9)')  X(1) 

WRITEl  ST0RE2.  »(f25.9)')  X(2) 
CMD»'run.ibhvg2  '//ST0RE1//'  '//ST0RE2 


c  UNIX  system  call 

°  CALL  SYS(CND) 

0PEN(9,f i le«' ibhvg2.out' ) 

REUIN0(9) 

REA0(9,*)  P 
READ(9,*)  V 
READ(9,*)  C 

write(*,'(f10.7,f10.7,f10.1,f10.1,f10.3)')  X<1),X(2),V,P,C 

V«-V 

RETURN 

END 
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SYS 


/*  this  C  function  is  used  to  call  a  UNIX  shell  conmand  (arbitrarily 
set  to  120  characters)  and  is  included  because  F77  does  not  have 
a  similar  facility 

*/ 


void  SYS(strng) 
char  *strng; 

^  *(strng>119)a'  •;  /*  force  the  string  to  be  null  terminated  V 


> 


sy8tem(strng); 


The  UNIX  shell  script,  run.ibhvg2,  executes  the  IBHVG2  code.  Prior  to  this,  it  modifies  an 
input  deck  called  ibhvglin  to  account  for  new  input  arguments  (/  and  D).  Muzzle  velocity, 
maximum  pressure,  and  charge  weight  are  passed  back  into  the  main  program  via  the  output 
file  ibhvglout.  This  is  rather  awkward,  but  the  method  is  very  general;  it  should  work  with 
most  any  operating  systems.  UNIX  provides  a  more  elegant  approach  through  the  use  of 
“pipes”  which  are  described  briefly  in  the  appendix. 


run.ibhvg2 


^  ----------aas=asssssarssssa==ss==ss========================= - — — — — — — — — — — 

#  Purpose:  Filter  for  IBHVG2  Input  and  Output: 

#  Filename:  run.ibhvg2 

#  Date:  22  July  1987 

#  sssaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa; 

# 

#  GENERAL  SCRIPT  INFORMATION: 

if  - 

#  Shell  script  input  argunents:  SI  -  DIAM  (  grain  diameter  in  main  propellant  deck) 

#  S2  -  LEN  (  grain  length  in  main  propellant  deck) 

# 

#  IBHVG2  input  filename:  ibhvg2.in 

# 

#  Shell  script  output  file:  ibhvg2.out 

#  -  ibhvg2.out  contains  the  values  for  maximum  pressure  (PMAX),  muzzle  velocity  (VMU2), 

#  and  the  charge  weight  of  the  main  propelling  charge  (CHWT),  respectively,  on 

#  separate  lines. 

# 

#  ALGORITHM: 


#  - 

# 

#  This  UNIX  shell  script  will  insert  the  values  for  LEN  and  DIAM  and  calculate  and  insert  the 

#  values  for  CHWT  and  web  (WEB)  into  the  IBHVG2  input  file.  This  input  file  is  directed 

#  into  the  IBHVG2  program.  The  values  for  VMUZ,  PMAX,  and  CHWT  are  filtered  from  output 

#  of  IBHVG2  and  placed  into  this  script's  output  file. 

# 

#  Charge  weight  will  be  determined  as  a  function  of  its  diameter,  length,  and,  in  some 

#  cases,  as  a  function  of  the  maximum  allowable  chamber  pressure,  using  the  following  algorithm: 

# 

#  1)  IBHVG2  will  adjust  CHWT  so  that  a  PMAX  of  79000  psi  is  obtained  (by  using 

#  the  PMAX  deck  option  in  IBHVG2). 

#  2)  The  maximum  loadable  charge  weight  (CHWT  MAX)  is  calculated  as  a  function  of  LEN 

#  and  DIAM  (and  system  constants)  from  an  empirical  formula  provided  by 

#  Kevin  Resnik. 

#  3)  If  the  CHWT  obtained  from  Step  1.  is  less  than  or  or  equal  to  CHWT^MAX, 

#  the  results  obtained  from  the  IBHVG2  run  in  Step  1.  are  used  as  the  output. 

# 

#  If  the  CHWT  obtained  from  Step  1.  is  greater  than  the  CHWT  MAX,  IBHVG2  will 

#  be  rerun  using  CHWT  MAX  as  the  actual  charge  weight  for  the  given  set  of  inputs. 
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NOTES  OF  CAUTION: 


This  shell  script  is  NOT  bullet-proof.  The  IBHVG2  input  file  must  be 
set  up  with  the  following  restrictions: 

a)  The  main  propellant  deck  must  be  placed  in  the  second  SPROP  deck. 

b)  Each  parameter  varied  by  this  shell  script  may  SOLELY  occupy 
one  line  of  the  1BHVG2  input  file. 

c)  For  each  parameter  varied  by  this  shell  script,  there  must  be 
only  one  space  before  and  after  the  equals  sign. 

Section  1:  Run  IBHVG2  for  for  CHUT  which  results  froa  a  PfMX  deck. 


cp  ibhvg2.in  inputi 


#  Split  input  file  into  sections  so  that  decks  can  be  easily  altered  using  sed  coomand. 

csplit  -k  -s  -f  part  inputi  '/\*PROP/'  <9>  2>/dev/null 

m  partOO  header .deck 
mv  partOI  propl.deck 
fflv  part02  tinp2.deck 
cat  part0[3-91  >  tail. deck 

#  Calculate  web  size  in  order  to  input  parameter  into  1BHVG2  input  deck.  (  PO  »  perf  diameter  ) 
P0*0.015 

WEB»'echo  "(SI  -  5.0*0.015)/6.0"  |  be  -I' 

#  Substitute  grain  parameters  into  main  propellant  deck. 

cat  tmp2.deck  I  sed  's/DIAh  »  .VOIAM  »  'S1'/p'  | 
sed  's/LEH  »  .VLEN  »  'S2'/p'| 
sed  's/UEB  »  ."/WEB  »  'SWEB'/p'  >  prop2.deek 

#  Combine  input  decks  into  one  IBHVG2  input  file. 

cat  header. deck  propl.deck  prop2.deck  tail. deck  >  inputi .deck 

#  Redirect  input  file  into  IBHVG2. 

IBNVG2  <  inputi .deck  >  output  2>/dev/null 

#  Filter  out  pertinent  data  from  output  file. 

PMAX»‘grep  "BREECH  PRESS"  output  |  awk  '{  printf  "X10.3f",  S4  >»' 

VHUZ>*grep  "VELOCITY"  output  |  awk  '{  printf  "X30.20f",  $4  >'» 

CHUT«*sed  -n  '/-  CHARGE  2  CHARGE  3  -/p'  output  1  grep  WEIGHT  |  awk  '{  print  S8  }'» 


Section  2:  Determination  of  maxinua  loadable  charge  weight  using  Kevin  Resnik's 
eapiricsl  charge  loading  formula  for  19  -  perf  cylindrical  propellant. 


DIAM.S1;  LEH-$2;  PO-0.015;  VC-SU.O;  RH00.06 
PERCEHT»*eeho  -13.125LEH  ♦  68.125  |  be  -I* 

RHOLOAO»*echo  "((SRHO)*(SOIAM*$OIAM  -  19.0*SPO*SPO))/(SOIAM"SOIAM)"  |  be  -I* 
CHWT  MAX-'echo  "SPERCENT*SRHOLOAD*SVC/100.0  "  |  be  -I' 
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»  -  -  -  — -  . 

#  Section  3:  If  CHUT  from  PMAX  deck  >  CHUT^MAX,  then  IBHVG2  is  rerun  using 

#  CHUT  MAX  as  the  actual  CHUT. 

# 

#  -  - — - - - 

#  Compare  CHUT  to  CHUT.MAX. 

if  C  'echo  “SCHUT"'  -gt  'echo  *'$CHUT_KAX**'  1 
then 

#  Prepare  new  IBHVG2  input  deck  with  CHUT^NAX  as  the  actual  charge  weight  used, 
sad  's/CHUT  «  .VCHWT  «  'SCHUT.MAXVp"  prop2.deck  >  prop2a.deck 

#  Remove  SPMAX  option  from  IBHVG2  input  deck. 

sed  'sASPMAXASCOMM/'  tail. deck  >  taila.deck 

#  Combine  input  decks  into  an  IBHVG2  input  file. 

cat  header. deck  propl.deck  prop2a.deck  taila.deck  >  input2.deck 

#  Redirect  input  file  into  IBHVG2  program. 

IBHVG2  <  input2.deck  >  output  2>/dev/null 

#  Filter  out  pertinent  data  from  IBHVG2  output  file. 

PMAX='grep  "BREECH  PRESS"  output  I  awk  •i  printf  "%10.3f",  S4  >'' 

VMUZ='grep  "VELOCITY"  output  |  awk  •i  printf  "%30.20f",  S4  >'' 

CHUT* 'sed  -n  CHARGE  2  CHARGE  3  -/p'  output  |  grep  UEIGHT  |  awk  'C  print  %S  >'' 

fi 

ih - Remove  Temporary  Files  - 

rm  *.deck  output*  part*  input* 

#  -  FINAL  OUTPUT  - 

#  write  output  data  to  ibhvg2.out 

echo  SPMAX  >ibhvg2.out 
echo  SVMUZ  »ibhvg2.out 
echo  SCHUT  »ibhvg2.out 

#  - 


The  script  modifies  the  input  file  shown  on  the  following  page  and  executes  IBHVG2  with 
these  incorporated  changes.  It  should  also  be  noted  that  the  feasibility  of  the  charge  weight  is 
examined  inside  the  script.  If  the  charge  weight  corresponding  to  the  maximum  design  pres¬ 
sure  exceeds  what  is  physically  possible  to  include  within  the  chamber  (see  Figure  2),  then  the 
maximum  loading  weight  will  be  used  instead. 


10 


T 


ibhvg2an 


SCOMM 

$HEAT 

SGUN 

$PROJ 

SRESI 

SINFO 

SRECO 

$PRIM 

SPROP 


SPROP 


INPUT  FILE  FOR  OPTIMIZATION  TEST  CASE 


TSHL 

TUAL 


0.00450 

293 


CSHL 

HO 


NAME  «  'GERMAN  120MM  GUN' 
LAND  ■  4.724  G/L  «  1. 


1848 

0.0648 

CHAM 

TRAV 


RSHL  a  0.284 
HL  «  1 


584 

187.11 


GRVE  a  4.724 
TWST  a  99 


NAME 

«  'SLUG' 

PRUT  a  18.00 

NPTS 

«  4 

AIR  a  1 

TRAV 

«  0, 

.8,  3.0,  187 

PRES 

«  100, 

2500,  100,  100 

RUN 

a  'OPTIMIZE'  CELT  «  5E-5 

DELP  a  5E-5 

EPS 

s  0.002 

GRAD 

a  2 

POPT  a  1,0, 1,0,0 

SOPT  a  0 

CONP 

*  0 

NAME 

a  'HOME' 

RECO  a  0 

RCUT  a  0 

$  Primer  Action  -  4.96  percent 

of  total  benite 

primer 

charge 

NAME 

«  'BENITE' 

CHUT  a  0.00347 

GAMA  a  1.25 

FORC 

«  212500 

COV  a  30 

TEMP  a  2000 

S  Ignitor  Action  -  95.04  percent 

of  total  benite 

primer 

charge 

NAME 

*  'BENITE' 

CHUT  a  0.06653 

GRAM  a  'CORD' 

IGNC 

s  0 

RHO 

s  0.06 

GAMA  a  1.25 

FORC  a  212500 

EROS 

3  0.000 

COV 

«  30 

TEMP  a  2000 

ALPH  a  0 

BETA 

»  27 

LEN 

»  9.998 

DIAM  a  0.078 

S  Main 

Propelling  Charge 

NAME 

«  'MAIN' 

GRAM  a  '19P' 

GAMA  a  1.257 

FORC 

*  39500( 

COV 

«  31.22 

TEMP  a  3049 

EROS  a  0.0000 

RHO  • 

>  0.06 

PO  a  0.015  NTBLa5 

PR4La  2000,  4000,  10000,  15000,  25000 

BR4La  0.292,  0.856,  1.904,  2.143,  3.394 

CHWT  a  19.10  %  ^  ^ 

LEN  ■  0.408  *  Remember:  One  and  only  one  space  before  and 

DIAM  ■  0.357  *  and  after  the  equals  sign  for  the 

WEB  «  0.047  t  variables  CHUT.LEN.OIAM.  &  WEB. 


JPROP  $  Combustible  Case  -  Modeled  as  a  deterred  propellant. 

NAME  »  'FMC  CASE'  CHUT  »  1.41  GRAM  «  '1PF'  RHO  » 
COV  »  27.927  TEMP  a  1610  EROS  »  0.00  GAMA  * 
FRCPa  ,150000,150000  FRCE*  ,150000,150000  FRCL<4)=200000 


0.04 

1.258 


IGNS(3)a2 
NTBLa2 

PR2P>1000, 10000 
PR2E>1000, 10000 
PR4La1 000. 10000 
LEM  -  18 


THRS(3)s200 
DEPP*  .  .015 
BR2P-.5,2.4 
BR2E«.5,2.4 
BR4La.5,10 
DIAM  ■  6.17 


.015  ,  .0155 


DEPEa 
.0155 

PR3P-1000, 10000  BR3Pa.5,2.4 
PR3Ea 1000, 10000  BR3Ea.5,2.4 


PO 


6.01 


UI 


.08 


SPROP 


S  Combustible  Adapter 


MAME  a  'KRAFT  CASE'  CHWT  a  .21 

COV  a  9.883  TEMP  a  1054 

LEM  a  3.4  DIAM  -  6.17 


GRAM 

EROS 

PO 


'1PF' 

0.00 

6.01 


SPMAX 


SEMD 


VARY  a  'CHWT'  MTH 

PMAX  a  79000  EPS 


a  2 
a  5 


TRY1  a  18.0 
LOOP  a  100 


RHO  a  0.04 
ALPH  a  1 
UI  a  .08 


TRY2  a  19.0 


GAMA 

BETA 


1.2734 

0.00001 


FORC  a  95726 
IGMC  a  0 
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Tables.  Optimization  Results 


Iteration 

Grain 

Diameter 

(inches) 

Grain 

Length 

(inches) 

Muzzle 

Velocity 

(fiA) 

Maximum 

Pressure 

(psi) 

Charge 

Weight 

(lbs) 

1 

03570000 

0.4080000 

5344.6 

78996.0 

19.891 

2 

0J681S63 

0.4080000 

5138.8 

71468.0 

19.842 

3 

03458437 

0.4080000 

5368.1 

79000.0 

19.423 

4 

03570000 

0.4207501 

5347.4 

im9.o 

19.930 

5 

03570000 

03952499 

5340.1 

79000.0 

19350 

6 

036S1S63 

0.4207501 

51073 

70137.0 

19.789 

7 

03681563 

03952499 

51713 

72873.0 

19.895 

8 

03458437 

03952499 

5363.8 

79000.0 

19385 

9 

03458437 

0.4207501 

5368.2 

79000.0 

19.460 

10 

03573486 

0.4080000 

53423 

78998.0 

19.905 

11 

03570000 

0.4083984 

5344.8 

789%.0 

19.892 

12 

03370761 

03953711 

5263.1 

71942.0 

18.945 

13 

03470381 

0.5016855 

5388.9 

79000.0 

19.702 

14 

03473TJ0 

03016855 

5ZS7J 

79000.0 

19.717 

15 

03470381 

03021755 

5389.0 

79000.0 

19.703 

16 

03349719 

0.6175594 

5256.1 

71512.0 

18.844 

17 

03410050 

0.5596225 

5406.1 

78996.0 

19.529 

18 

03413380 

0.5596225 

5404.1 

78995.0 

19.545 

19 

03410050 

0.5601690 

5406.2 

789%.0 

19.530 

20 

03361544 

0.6030530 

5261.0 

71913.0 

18.909 

21 

03385797 

03813377 

5411.8 

79003.0 

19.447 

22 

03389104 

0.5813377 

5412.8 

78999.0 

19.464 

23 

03385797 

0.5819055 

5411.9 

79004.0 

19.448 

24 

03390746 

0.5856780 

5248.6 

71472.0 

18.993 

25 

03388272 

0.5835079 

5412.8 

78999.0 

19.463 

26 

03391580 

03835079 

5412.0 

79000.0 

19.479 

27 

03388272 

0.5840777 

5254.7 

71753.0 

18.999 

28 

03388284 

03835007 

5412.8 

78999.0 

19.463 

29 

03388296 

03834935 

5412.8 

78999.0 

19.463 

30 

03388320 

03834790 

5412.8 

78999.0 

19.463 

31 

03388368 

03834502 

5412.9 

79000.0 

19.463 

32 

03388464 

0.5833925 

5412.9 

79000.0 

19.464 

33 

03391773 

03833925 

5412.0 

78998.0 

19.479 

34 

03388464 

0383%23 

5254.6 

71751.0 

19.000 

35 

03388465 

03833916 

541Z9 

79000.0 

19.464 

36 

03388466 

03833908 

5412.9 

79000.0 

19.464 

37 

03391774 

03833916 

5412.0 

78997.0 

19.479 

38 

03388465 

03839614 

5254.6 

71751.0 

19.000 

39 

03388465 

03833916 

541Z9 

79000.0 

19.464 

40 

03391774 

03833916 

5412.0 

78997.0 

19.479 

41 

03385156 

03833916 

5411.6 

79001.0 

19.441 

42 

03388465 

03839614 

5254.6 

71751.0 

19.000 

43 

03388465 

03828219 

5412.8 

79000.0 

19.463 

44 

03389031 

03832820 

5413.0 

78997.0 

19.466 

45 

03389598 

03831724 

54133 

78998.0 

19.469 

46 

03390731 

0.5829533 

5411.8 

78999.0 

19.474 

47 

03392908 

0.5831724 

5412.1 

78995.0 

19.484 

48 

03386288 

03831724 

5412.1 

79003.0 

19.452 

49 

03389598 

0.5837420 

54133 

78998.0 

19.470 
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The  results  of  the  optimization  are  shown  in  Table  3.  As  can  be  seen,  the  program  finds  the 
crest  region  and  then  searches  slowly  due  to  the  flat  slope  near  the  top.  Clearly,  this  would  be 
very  tedious  to  perform  by  hand.  Increasing  the  number  of  variables  to  be  optimized  would 
quickly  render  the  problem  intractable  without  a  computational  scheme  and  a  fast  computer. 
To  provide  some  feel  for  the  computational  effort,  this  run  requires  about  five  minutes  on  a 
supercomputer  and  about  two  hours  on  a  minicomputer. 


B.  Two  Thermodynamic  Examples 

The  BLAKE  code  computes  thermodynamic  and  equilibrium  information  for  hot  propelling 
gases  in  a  ballistic  environment.  There  are  many  circumstances  where  running  the  BLAKE  code 
in  an  iterative  fashion  would  be  useful  Our  first  case  modifies  the  concentration  of  a  two- 
component  propellant  to  search  for  the  combination  of  maximum  total  impetus. 


Table  4.  Program  Organization  for  Case  1 


Figure  1 

Example  1 

Description 

MAIN 

MAXIMP 

main  program  for  optimization 

SUBl 

FI 

user-written  routine  called  by 
ZXLSF  (IMSL  routine) 

SYS 

sys 

C-code  to  perform  system  call 

Commands 

nm.blake 

operating  system  commands  to 

execute  “black  box” 

MAIN2 

BLAKE 

“black  box”  code 

The  idea  is  virtually  the  same  as  before.  The  main  program  sets  up  the  problem,  an  IMSL  rou¬ 
tine  does  the  work  which  executes  a  shell  script.  The  IMSL  routine  for  this  problem,  ZXLSF, 
performs  a  one-dimensional  minimization  of  a  smooth  function  using  a  safeguarded  quadratic 
interpolation.  The  results  are  shown  graphically  in  Figure  4.  The  coding  follows. 


MAXIMP 


program  maximp 


external  f 

integer  ier,  mexfn 

real  f,  xacc,  x,  step,  bound 


X  >  50.0 
step  s  *40.0 
bound  «  49.0 
xacc  «  0.01 

maxfn  «  50 


#  first  guess,  10.0*50^step 

#  1.0  <■  X  <■  99.0 

#  accuracy 

#  maxinui  function  evaluations 


call  zxlsf(f1,  X,  step,  bound,  xacc,  maxfn,  ier) 
write  (6,  '<a,f11.5,a,i3)^  )'X«',x,MER«',ier 
end 


Fl 


real  function  f1(x) 
c 

c  call  a  Unix  shell  script  which  executes  the  BLAKE 

c  computer  code  with  input  parameter  x  (X  concentration  of 

c  polyethylene)  and  returns  the  value  impet  (propellant  impetus) 
c 


real  x,  polyet,  impet 
character*60  cmd 
character*20  storel 
c 

c  using  internal  writes  to  convert  the  floating  point 
c  numbers  to  a  character  string 
c 

writeC  storel,  '(g15.6)0  x 
cmdB'run.blake  V/storel 
c 

c  Unix  system  call 

c 

call  SYS(cmd) 

open(9,f i les'blake.out' ) 
rewind(9) 

read(9,*)  polyet, impet 

write(*,'(a,f10.2)^)  '  concentration  of  poly  eth.*:  ^polyet 

write<*,'(a,f10.1//)^)  •  impetus  is:  ^impet 

f*-impet 

return 

end 


blake.in 


PRL,CON,0,PAG,0 

TITLE,  POLYETHYLENE  PLUS  70X  H202 
FOR,  POLYET,  -UOOO,  C,2,H,4 
FOR,OXID70,-1558.28E6,0,43705,H, 56295 
COM,  POLYET,  8.0000  ,  0X1070,  92 
UNI,  ENG 
GUN,  0.2,10,0 
STOP 
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nin.blake 


«  BLAKE  Driver: 

P0LY»$1  #  Oxidant  is  (  100X  -  Polyethylene  concentration  ) 

OXY«*  echo  "100.0  -  BPOLY  "  |  be  * 

sed  'S/COM,  POLYET,  8.0000/COM,  POL YET,  'SPOLY'/p'  <blake.in  | 
sed  's/OXIDTO,  92/OXID70,  »SOXY'/p'  >teflp 

BLAKE  <  tenp  >  output  2>/dev/null 


#  inpetus  is  the  2nd  from  the  last  line,  in  colunn  5 
IMPETUS«’tail  -2  output  |  sed  -n  Ip  |  auk  '€  printf  "*10.3f",  $5  > 


# 

#  -  FINAL  OUTPUT 

# 

echo  SPOLY  >blalceeOUt 
echo  $IMPETUS  »blake.out 

# - 


Figure  4,  Optimization  of  Impetus 


The  second  case  again  adjusts  the  combination  but  now  we  are  looking  for  the  lowest  con¬ 
centration  of  polyethylene  which  produces  no  solid  carbon  as  a  reaction  product  (fouls  gun 
mechanisms).  The  search  algorithm  is  a  simple  bisection  method  called  directly  from  the  main 
program. 


Table  5.  Program  Organization  for  Case  2 


Figure  1 

Example  1 

Description 

MAIN 

BRACKET 

main  program  for  optimization 

SUBl 

F2 

user-written  routine 

SYS 

sys 

C-code  to  perform  system  call 

Commands 

run.blake.2 

operating  system  commands  to 

execute  “black  box” 

MAIN2 

BLAKE 

“black  box”  code 

BRACKET 


c  main  program 
external  f2 

real  f2,  xO,  x1,  x,  tol,  fx 

tol 30.005 

x0«1.0 

x1»99.0 

c  repeat  until  differences  are  <  tol 
10  continue 

xs0.5*(xUx0) 

fx*f(x) 

if  <  fx  .It.  tol  )  then 
xO^x 
else 
xlsx 
endif 

if  (  (xl-xO)  .gt.  tol  )  go  to  10 

urite  (*,  ^(a.f6.2.a,f15.5)'  )  'X  *  \  x/  C(s)  *',fx 
end 
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F2 


real  function  f2(x) 

c  call  a  Unix  shell  script  which  execute  the  blake 
c  computer  code  with  input  parameter  x  (X  concentration  of 
c  polyethylene  and  returns  the  value  CS  (solid  carbon) 
c 

real  x.  polyet,  cs 
character*60  cmd 
characterize  storel 

c  using  internal  writes  to  convert  the  floating  point 
c  numbers  to  a  character  string 
c 

write(  storel,  '(g15.6)')  x 
cmd* 'run. blake. 2  '//storel 
c 

c  Unix  system  call 

call  SYS(cmd) 

open(9,file«'blake.out') 

rewind(9) 

read(9,i)  polyet.cs  ^  ,14. 

write(i,'(a,f10.3)')  '  concentration  of  poly  eth.*:  ',polyet 

write(i,'(a,f10.3//)')  '  solid  carbon  is:  ',cs 

f=cs 

return 

end 


run.blake.2 


#  BLAKE  Driver: 

P0LY«*1  #  Oxidwit  is  (  100X  -  Polyethylene  concentration  ) 

OXY.‘  echo  "100.0  -  SPOLY  «  |  be  ' 

sed  'S/COM,  POLYET,  8.0000/COM,  POLYET,  'SPOLY'/p'  <blake2.in  | 
sed  's/OXIDTO,  92/0X1070,  '$OXY'/p'  >te)rp 

BLAKE  <  temp  >  output  2>/dev/null 

Cs»*  grep  'C(S)  SOLID'  output  |  awk  '<  printf  "X15.5f",  $4  >'* 

echo  SPOLY  >blake.out 

echo  SCs  »blake.out 


blake2.in 


TITLE,  POLYETHYLENE  PLUS  70*  H202 

FOR,  POLYET,  -14000,  C,2,H,4 

FOR,OXID70,-1558.28E6,0,43705,H,56295 

COM,  POLYET,  8.0000  ,  OXID70,  92 

UNI,  ENG 

GUN,  0.2,10,0 

STOP 
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Figure  5.  Maximum  Concentration  without  Solid  Carbon 


in.  FURTHER  EXTENSIONS 


It  is  hoped  that  cases  described  in  the  previous  section  are  sufficient  to  illustrate  details 
of  the  method  yet  simple  enough  for  the  reader  to  comprehend  each  problem  without  undue 
detail.  Practical  problems  would  have  main  programs  of  several  hundred  lines  and  numerous 
rails  to  various  subroutine  libraries.  This  section  is  included  to  briefly  discuss  more  elaborate 
applications  of  this  technique  which  may  highlight  its  strengths  and  versatility. 

A.  Parallel  Processing 

Many  computational  schemes  require  a  number  of  similar  calculations  which  are  uncoupled 
from  each  other.  As  an  illustration,  ab  initio  chemistry  problems  typically  require  large 
numbers  of  individual  integrals  to  be  evaluated.  Running  these  independent  sub-problems  on 
multiple  processing  elements  in  parallel  can  be  quite  attractive.  The  taxonomy  of  parallel  com¬ 
putation  is  best  left  to  the  literature  [8];  however,  the  different  approaches  are  characterized  by 
1)  the  number  of  independent  paths  to  memory  with  data  streams  which  can  be  shunted  in 
parallel,  and  2)  the  number  of  independent  instructions  which  can  be  executed  in  parallel.  The 
multiple  instruction  stream,  multiple  data  stream  (MIMD)  machine  is  rapidly  becoming  the 
most  cost-effective  design  because  the  processors  tend  to  operate  independently.  This  permits 


18 


an  architecture  of  inotpensive,  off-the-shelf  computers.  The  tradeoff  is  usually  a  more  complex 
algorithm  in  order  to  effectively  harness  the  multiple  processors.  The  problems  which  work 
best  are  those  which  can  be  decomposed  into  distinct  processes  conditioned  to  run  on  separate 
computers  simultaneously.  That  is,  they  should  require  no  communication  at  intermediate 
steps.  All  the  better  if  the  processes  execute  in  a  predictable  manner.  This  alleviates  the  over¬ 
head  of  synchronizingthe  results  and  efficiendy  loading  the  multiple  processors. 

The  type  of  problems  discussed  in  this  paper  are  ideal  for  this  situation.  Separate  “black 
box”  runs  could  be  pawned  simultaneously  and  processed  independently  via  distinct  comput¬ 
ers.  The  optimization  example  with  IBHVG2  could  easily  exploit  this  strategy.  A  Hessian 
matrix  is  computed  which  requires  several  predetermined  IBHVG2  runs  to  be  executed.  A 
serial  machine  would  do  one  after  the  other,  while  a  parallel  configuration  could  do  all  the  cal¬ 
culations  for  the  matrix  at  the  same  time.  This  is  illustrated  in  Figure  6.  Furthermore,  the  exe¬ 
cution  times  of  each  run  of  IBHVG2  in  this  circumstance  will  be  nearly  identical.  Although  the 
complexity  of  the  main  program  would  increase,  the  burden  should  not  be  too  oppressive  and 
the  payoff  could  be  substantial 
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Figure  6.  Potential  Speedup  due  to  Parallel  Processing. 
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B.  Expert  Systems 

The  terms  artificial  intelligence  and  expert  systems  have  become  popular  buzz  words  of  the 
1980’s.  Most  reports  and  proposals  which  use  them  rely  more  on  flash  than  substance.  This  is 
particulary  true  when  the  applications  wander  outside  the  realm  of  computer  science.  The 
hopes  for  expert  systems  are  usually  the  modem  day  version  of  a  “free  lunch” — computers  with 
some  magical  software  are  suppose  to  solve  formidable  problems  without  the  need  of  human 
experts.  The  question  is:  who  writes  the  software?  No  successful  e^ert  system  has  avoided 
the  active  involvement  of  the  experts.  This  is  why  there  are  so  few  expert  systems  which 
accomplish  significant  tasks.  By  definition,  experts  are  highly  skilled  individuals  in  their  chosen 
disciplines.  Convincing  them  to  allocate  sizable  chunks  of  time  to  produce  a  complex  computer 
program  can  be  both  difficult  and  very  expensive.  Most  commercial  ®qpert  systems  are  simply 
elaborate  database  programs  with  little  deductive  reasoning  capability. 

This  may  well  change  in  the  next  decade  or  two.  A  recent  resurgence  in  artificial  intelli¬ 
gence  has  been  fostered  by  the  dramatic  increase  in  computing  resources  and  some  novel 
hardware  concepts.  Neural  networks  and  massively  distributed  computing  are  two  of  the  latest 
topics  under  scrutiny.  Replacing  or  augmenting  experts  with  software  is  certainly  a  long-range, 
but  laudable  goal.  The  commercial  impact  of  synthesizing  human  reasoning  and  experience 
gathering  would  be  enormous.  Although  researchers  clearly  have  many  obstacles  to  overcome, 
the  potential  payoff  is  sufficient  to  sustain  interest  for  many  years  to  come. 

In  the  meantime,  we  suggest  a  less  ambitious  approach.  Rather  than  trying  to  replace  the 
expert,  it  is  probably  more  fruitful  to  merely  automate  routine  manipulations.  It  is  possible  to 
utilize  the  working  computational  codes  of  the  experts  and  couple  them  to  another  code  which 
could  simulate  some  simple  thinking  processes.  We  prefer  to  call  this  portion  “automated  rea¬ 
soning”  to  perhaps  convey  the  more  limited  scope.  There  are  several  realizable  benefits  in  a 
design  situation.  Reduced  labor  costs  are  possible  by  eliminating  most  of  the  user  interaction 
during  the  design  iteration  cycles.  Enhanced  designs  are  possible  because  the  computer  can  sys¬ 
tematically  examine  numerous  designs  far  beyond  the  patience  and  endurance  of  a  human 
operator.  An  intelligent  interface  can  minimize  errors  by  screening  information  for  consistency, 
satisfy  constraints,  etc.  If  nothing  else,  the  computational  modules  can  be  made  less  cumber¬ 
some  by  driving  them  with  a  user-friendly  environment.  A  simple  mock-up  is  illustrated  in  Fig¬ 
ure?. 

There  is  a  fundamental  difference  between  this  use  of  artificial  intelligence  and  the  applica¬ 
tions  one  sees  in  magazines  or  on  television.  The  popularized  programs  are  knowledge-based  in 
that  they  contain  a  very  complex  database  of  knowledge  distilled  from  experts.  Elaborate  and 
sophisticated  logical  inferences  require  fast  manipulation  of  the  database  which  may  be  filled 
with  qualitative  and  subjective  information.  An  example  of  this  type  program  is  an  expert  sys¬ 
tem  for  medical  diagnosis.  Our  approach  is  compute  intensive  as  opposed  to  logic  intensive. 
The  “intelligence”  resides  mainly  in  the  experts’  computational  codes.  The  computer  runs  them 
much  like  the  human  would — just  thousands  of  times  faster.  It  is  analogous  to  observing  an 
experiment  for  quantitative  information.  In  fact,  the  computer  models  usually  are  a  distillation 
of  experiments  themselves  and  they  simulate  certain  features  which  we  understand.  As  opposed 
to  knowledge-based,  we  refer  to  this  type  of  expert  system  as  artafysis-based. 
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Figure  7.  Proposed  System  for  Charge  Design 


21 


Optimization 

Analysis 


Expert 

Systems 


Code  Se 

lection 

_ i 

i 

^Reconfigure  cbargr^ ; 
.^based  on  logic| 


. t 


Lumped'parameter  Codes 
•BLAKE 
•  IBHVG2 


No 


Distributed-parameter  Codes 

•  NOVA 
•TDNOVA 

•  XKTC 

- f - 


FFT 


Figure  8.  Sample  Design  Strategy 


22 


Logical  inferences  could  be  made  without  introducing  the  specter  of  artifical  intelligence 
(AI).  The  reasoning  section  could  be  implemented  in  a  FORTRAN  code  just  by  constructing 
intricate,  nested  If-Then-Else  statements.  However,  AI  languages  such  as  LISP  and  PROLOG 
are  far  more  suitable  to  handle  logic.  The  key  features  which  make  these  languages  ideal  is  the 
ability  to  manipulate  program  source  code  as  data  (self-modifying),  the  capability  to  manipulate 
symbols  as  well  as  numbers,  and  finally,  (fynamic  allocation  of  storage.  Together,  these  advan¬ 
tages  produce  concise,  straight-forward  programs  which  can  be  molded  to  applications  easily. 
Most  expert  systems  are  rule-based.  This  means  that  the  execution  and  communication  of  dif¬ 
ferent  modules  can  string  together  much  like  theorems  in  mathematics.  Furthermore,  there  is 
usually  a  way  of  tracing  through  the  series  of  rules  to  examine  how  the  program  arrived  at  a 
particular  result.  This  is  called  back-tracking.  It  is  extremely  useful  for  debugging  and  verifying 
the  logic.  We  will  not  go  into  the  details  of  AI  logic  programming  here.  The  main  point  is  that 
coupling  an  AI  language  with  FORTRAN  analysis  codes  can  be  effectively  implemented  using 
the  techniques  described  in  the  report. 

Let  us  conjure  up  a  possible  ballistics  example.  Suppose  one  needed  to  design  a  propulsion 
system  for  a  nuclear  round.  Besides  the  standard  performance  requirements  of  zone  charges 
and  muzzle  velocities  with  range  overlap,  great  care  is  to  be  paid  to  pressure  waves  due  to  the 
nature  of  the  payload.  The  design  strategy  might  look  something  like  Figure  8. 

An  initial  charge  configuration  is  provided  as  input  to  both  BLAKE  and  IBHVG2  to  com¬ 
pute  and  optimize  the  gross  ballistics  of  the  round  (mean  pressure  and  velocity  history).  These 
codes  would  not  be  capable  of  predicting  pressure  waves,  however.  At  an  appropriate  point,  a 
distributed-parameter  code  such  as  NOVA  would  be  needed.  The  resulting  wave  behavior 
could  be  evaluated  by  computing  the  Fourier  transform  of  the  breech  pressure  history.  A 
power  spectrum  beyond  a  specified  norm  could  flag  an  unacceptable  design.  Depending  on  the 
results,  the  charge  could  be  modified  by  a  rule-based  expert  system.  For  example,  pressure 
waves  might  be  ameliorated  by  increasing  the  ullage  or  by  altering  the  ignition  system.  Radical 
changes  such  as  propellant  reformulation  would  require  going  all  the  way  back  to  the  lumped- 
parameter  codes,  while  perturbations  could  iterate  using  the  distributed  codes. 


IV.  CONCLUSIONS 


A  simple  method  of  tying  together  existing  computer  programs  has  been  both  demon¬ 
strated  and  explained.  The  implementation  involves  system  calls  to  the  operating  system. 
Although  the  examples  in  this  paper  have  focussed  on  the  UNIX  operating  system,  the  tech¬ 
nique  is  general  and  should  be  applicable  to  many  other  operating  systems.  Some  advantages 
to  this  approach  include  Increased  modularity,  added  productivity,  and  the  potential  for  more 
elaborate  applications.  It  is  hoped  the  three  examples  provided  in  the  text  have  demonstrated 
the  first  two  virtues.  The  last  advantage  was  only  discussed  and  is  still  to  be  proven.  A  convinc¬ 
ing  application  will  be  the  subject  of  a  future  technical  report. 


23 


REFERENCES 


[1]  R.  Anderson  and  K.  Fickie,  “IBHyG2— A  User’s  Guide,”  BRL-TR-2829,  Ballistic 
Research  Lalx)ratory,Aberdeen  Proving  Ground,  MD,  1987. 

[2]  E.  Freedman,  “BLAKE— A  ThermotfynamicsCode  Based  on  TIGER:  User’s  Guide  and 
Manual,”  ARBRL-TR-0241 1,  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
MD,  1982. 

131  IMSL  Library  Reference  Manual,  Edition  9.2,  (IMSL  Inc.,  7500  Bellaire  Boulevard, 
Houston  TX  77036)  Volume  4,  Chapter  Z,  1984. 

[4]  G.  Strang,  Introduction  to  Applied  Mathematics,  Wellesl^-Cambridge  Press  (Welles¬ 
ley,  1986); 

[5]  W.H.  Press,  B.P.  Flannery,  SA.  Teukolsky,  W.T.  Vetterling,  Numerical  Recipes,  Cam¬ 
bridge  University  Press  (New  York,  1986). 

161  P.E.  Gill,  W.  Murray,  M.  H.  Wright,  Practical  Optimization,  Academic  Press  (New 
York,  1981). 

(71  E.  Polak,  Computational  Methods  in  Optimization,  Academic  Press  (New  York,  1971). 

(8]  G.  Rodrigue,  Parallel  Computations,  Academic  Press  (New  York,  1982). 


24 


APPENDIX 


As  briefly  mentioned  in  the  report,  the  UNIX  operating  system  provides  additional  ways  of  exe¬ 
cuting  system  commands  inside  an  active  process  (Le.,  a  separate  executing  computer  program). 
One  technique  would  be  to  use  pqjes.  Pipes  allow  transfer  of  data  streams  between  processes 
and  can  be  used  for  process  synchronization.  Their  most  useful  characteristic  for  inter-process 
communication  is;  streams  need  no  information  about  what  processes  are  at  the  other  end  of 
the  pipeline.  There  are  two  types  of  pipes.  The  simplest  is  a  system  pipe  shown  in  Figure  9 
(some  people  call  these  unnamed  pipes).  Basically,  they  work  much  like  files,  but  the  UNIX 
kernel  uses  direct  blocks  of  the  inode  for  greater  efficiency.  The  other  type  is  a  process  pipe 
(sometimes  called  a  named  pipe).  Figure  10  shows  typical  code  fragments  using  process  pipes. 

It  would  be  inappropriate  to  fully  discuss  the  mechanisms  and  added  vutues  of  using  pipes 
for  the  applications  described  in  this  paper.  The  discussion  would  be  rather  involved  and  to 
fully  comprehend  the  arguments  requires  some  background  in  UNIX  system  programming.  We 
mention  them  because  a  subroutine  call  to  a  “black  box”  process  using  pipes  could  more  closely 
emulate  a  normal  subroutine  call  written  in  a  standard  language  such  as  FORTRAN.  Pipes  can 
both  send  values  and  return  values  transparently  without  resorting  to  external  files.  They  would 
provide  an  elegant  interface  between  processes,  and  potentially  more  sophisticated  and  capable. 
We  relegate  this  idea  to  an  appendix  because  the  method  is  very  specific  to  UNIX.  See  your 
system  programmer  for  more  details. 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Tour  coonaents/answers  to  the  iteas/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received  _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 
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City,  State,  Zip 
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New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 
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City,  State,  Zip 
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